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I review the Higgs sector of the U{1)b-l extension of the minimal supersymmetric standard model 
(MSSM). I will show that the gauge kinetic mixing plays a crucial role in the Higgs phenomenology. 
Two light bosons are present, a MSSM-like one and a. B — L-like one, that mix at one loop solely 
due to the gauge mixing. After briefly looking at constraints from flavour observables, new decay 
channels involving right-handed (s)neutrinos are presented. Finally, it will be reviewed how model 
features pertaining to the gauge extension affect the model phenomenology, concerning the existence 
of R-Parity-conserving minima at loop level and the Higgs-to-diphoton coupling. 


INTRODUCTION 


The recently discovered Higgs boson is considered as 
the last missing piece of the standard model (SM) of par¬ 
ticle physics. Nonetheless, several firm observations uni- 
vocally call for its extension. Mainly but not limited to, 
the presence of dark matter, the neutrino masses and 
mixing pattern, the stability of the SM vacuum, the hi¬ 
erarchy problem. Supersymmetry (SUSY) has long been 
considered as the most appealing framework to extend 
the SM. Its minimal realisations (MSSM and its con¬ 
strained versions start however to feel considerable 
pressure to accommodate the recent hndings, especially 
the measured Higgs mass of 125 GeV. Despite not in 
open contrast with the MSSM, the degree of fine tun¬ 
ing required to achieve it is more and more felt as un¬ 
natural. In order to alleviate this tension, non-minimal 
SUSY realisations can be considered. One can either ex¬ 
tend the MSSM by the inclusion of extra singlets {e.g. 
NMSSM i) or by extending its gauge group. Concern¬ 
ing the latter, one of the simplest possibilities is to add an 
additional Abelian gauge group. I will focus here on the 
presence of an [/(1)b_l group which can be a result of an 
Eg, X Eg heterotic string theory (and hence M-theory) 

[^. This model, the minimal i?-parity-conserving B — L 
supersymmetric standard model (BLSSM in short), was 
proposed in [^, Q and neutrino masses are obtained via a 
type I seesaw mechanism. Furthermore, it could help to 
understand the origin of i?-parity and its possible spon¬ 
taneous violation in supersymmetric models as well 
as the mechanism of leptogenesis 

It was early pointed out that the presence of two 
Abelian gauge groups in this model gives rise to kinetic 
mixing terms of the form 


a^b 


( 1 ) 


that are allowed by gauge and Lorentz invariance [ll| , as 


^ For a review, see Ref. Q. 


pa,uv pb,ni' g^]-g gauge-invariant quantities by them¬ 
selves, see e.g. 12|. Even if these terms are absent at 
tree level at a particular scale, they will in general be 
generated by RGE effects ISj, [l^ . These terms can have 
a sizable effect on the mass spectrum of this model, as 
studied in detail in Ref. 1^, and on the dark matter, 
where several scenarios would not work if it is neglected, 
as thoroughly investigated in Ref. [l^. In this work, 
I will review the properties of the Higgs sector of the 
model. Two light states exist, a MSSM-like boson and a 
B — L-like boson. After reviewing the model, I will show 
that a large portion of parameter space exists where the 
SM-like Higgs boson has a mass compatible with its mea¬ 
sure, both in a “normal” (M//^ > Mh^ = 125 GeV) and 
in a “inverted” hierarchy {Mh^ < Mh^ = 125 GeV), also 
in agreement with bounds from low energy observables 
and dark matter relic abundance. The phenomenological 
properties of the two lightest Higgs bosons will be system¬ 
atically investigated, where once again the gauge mixing 
is shown to be fundamental. The presence of extra D- 
terms arising from the new U (1) b-l sector, as compared 
to models based on the SM gauge symmetry, has a large 
impact on the model phenomenology. They affect both 
the vacuum structure of the model and the Higgs sector, 
in particular enhancing the Higgs-to-diphoton coupling. 
Both these issues will be reviewed here, despite the latter 
is disfavoured by recent data 17| , to show model features 
beyond the MSSM. 


THE MODEL 


Eor a detailed discussion of the masses of all particles 
as well as of the corresponding one-loop corrections we 
refer to [l^ . Attention will be payed on the main aspects 
of the U(l) kinetic mixing since it has important conse¬ 
quence for the scalar sector. Eor the numerical investi- 
gat ions that will be shown, we used the SPheno version 
3, 1^ created with SARAH 23-24| for the BLSSM. For 
the standardised model definitions, see Ref. [2^, while 
for a review of the model implementation in SARAH, see 
Ref. [2^. This spectrum calculator performs a two-loop 
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TABLE I: Chiral superfields and their quantum numbers un¬ 
der Gsm < 8 ) U{1)b-l, where Gsm = {U{1)y < 8 ) SU{2)l (8 
SU{3)c) ■ 


RGE evaluation and calculates the mass spectrum at one 
loop. In addition, it calculates the decay widths and 
branching ratios (BRs) of all SUSY and Higgs particles as 
well as low-energy observables like {g — 2)^. We will dis¬ 
cuss the most constrained scenario with a universal scalar 
mass mo, a universal gaugino mass M 1/2 and trilinear 
soft-breaking couplings proportional to the superpoten¬ 
tial coupling (Ti = AoYi) at the GUT scale. Other input 
parameters are tan/3, tan/3', Mz', Tx, and Y^. They 
will be defined in the following section. The numerical 
study here presented has been performed by randomly 
scanning over the independent input parameters above 
described via the SSP toolbox [^, while low energy ob¬ 
servables such as BR(/i —>• ey) and BR(/r —3e) have 
been evaluated with the FlavourKit package [ 2 ^. Fur¬ 
thermore, during the scans all points have been checked 


with HiggsBounds-4.1.1 |29l - t32| . both in the “normal” 


hierarchy and in the “inverted" hierarchy case. 


Particle content and superpotential 


The model consists of three generations of matter par¬ 
ticles including right-handed neutrinos which can, for 
example, be embedded in 50(10) 16-plets. Moreover, 
below the GUT scale the usual MSSM Higgs doublets 
are present as well as two fields rj and fj responsible 
for the breaking of the U{1)b-l- The g field is also 
responsible for generating a Majorana mass term for 
the right-handed neutrinos and thus we interpret its 
B — L charge as its lepton number. Likewise is for 
7 /, and we call these fields bileptons since they carry 
twice the lepton number of (anti-)neutrinos. The quan¬ 
tum numbers of the chiral superfields with respect to 
U{1)y X SU{2)l X SU{3)c x U{1)b-l are summarised 
in Table m 


The superpotential is given by 

W = y:^ u- Qj Hu - Y^/ dt Qj Hd - Yi^ el L, 

+fj, Hu Hd + UJ-' Lj Hu - g fjg + Yl;P vl fj 

( 2 ) 

and we have the additional soft SUSY-breaking terms: 
■^SB = dY'MSSM — A^Ajj,/Ubb' — 

+ Tl^HuutL, + (3) 

i, j are generation indices. Without loss of generality one 
can take and to be real. The extended gauge 
group breaks to SU{‘i)c®U{l)em as the Higgs fields and 
bileptons receive vacuum expectation values (vevs): 

Hd= ^ (CTd + Vd+ i(j)d ), H°=^[au + VuY i(j)u) 

g= -^{ar^+Vr^+icjir^) , g = ^ [af) + Vff + icjyf)) 

We define tan/3' = in analogy to the ratio of the 

MSSM vevs (tan/3 = Vu/vd)- 


Gauge kinetic mixing 


As already mentioned in the introduction, the pres¬ 
ence of two Abelian gauge groups in combination with 
the given particle content gives rise to a new effect ab¬ 
sent in any model with just one Abelian gauge group: 
gauge kinetic mixing. This can be seen most easily by 
inspecting the matrix of the anomalous dimension, which 
for our model at one loop reads 


^ 167r2 




( 4 ) 


with typical GUT normalisation of the two Abelian gaMe 
groups, i.e. 1 / 3/6 for U{1)y and \/3j2 for U{1)b-l [3- 
Therefore, even if at the GUT scale the 17(1) kinetic mix¬ 
ing terms are zero, they are induced via RGE evaluation 
at lower scales. It turns out that it is more convenient 
to work with non-canonical covariant derivatives rather 
than with off-diagonal field-strength tensors as in eq. ©■ 


However, both approaches are equivalent 33|. Therefore 


in the following, we consider covariant derivatives of the 
form = d^ — iQ^GA where is a vector containing 
the charges of the field (j) with respect to the two Abelian 
gauge groups, G is the gauge coupling matrix 


G = \ 

\ gBY gBB 

and A contains the gauge bosons A = (A/, A^)^. 
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As long as the two Abelian gauge groups are unbroken, 
we have still the freedom to perform a change of basis by 
means of a suitable rotation. A convenient choice is the 
basis where qby = 0, since in this case only the Higgs 
doublets contribute to the gauge boson mass matrix of 
the SU{2)l 0 U{1)y sector, while the impact of t] and fj 
is only in the off-diagonal elements. Therefore we choose 
the following basis at the electroweak scale : 


is imposed, be., = ggL^^'^ and = 

= 0, at SUSY scale we get [l^ 

gBL = 0.548, (10) 

g ~ -0.147. (11) 


9yy = 

QYYQBB —QYBQBY 

= 9i 

y/ 9 ‘bb'^9\y 

q'bb = 

V 9b B + 9 BY — 

9bl 

9yb = 

QY BUBB-\-9BY QYY 

= 9 

y/ 9 bb~‘^9%y 

9'by = 

0 



( 6 ) 

(7) 

( 8 ) 

(9) 


When unification at some large scale (^ 2 • 10^® GeV) 


Tadpole equations 

The minimisation of the scalar potential is here de¬ 
scribed in the so-called tadpole method. We can solve 
the tree-level tadpole equations arising from the mini¬ 
mum conditions of the vacuum with respect to g,, g' 
and B^i. Using vl = v‘^ + n? and = v‘^ + we obtain 


ImP = i((25<7BL^^xCOs(2 ^')sec(2/3)-4(m |^^- [ol + f + (12) 

Bf, = -l(^ - 2ggg^vlcos{2l3') + + (^gf + f + gfjv^ cos{2l3fj tan(2/3) (13) 

\g\^ = l(^-2(^g%^vl+ml + m^^ + (2ml-2ml+ggs^v^cos{2l3)^sec{2l3')^ (14) 

Bf,' = I ( - 2gl^vl cos(2/3') -h 2ml - 2m? -h cos( 2/3)) tan(2/3') (15) 


Mz' — gBh'^x and, thus, we find an approximate relation 
between Mz' and g' 

Mz'^ ^ -2|/rf 

4(m? — m^ tan^ j3') — v'^gg^^ cos /3(1 -I- tan /?') 
2(tan2 (3'-1) 

(16) 

For the numerical results, the one-loop corrected equa¬ 
tions are used, which lead to a shift of the solutions in 
eqs. (fT^-(fT5]). 


The scalar sector 

In this model, 2 MSSM complex doublets and 2 bilep¬ 
ton complex singlets are present, yielding 4 CP-even, 2 
CP-odd, and 2 charged physical scalars. 

Concerning the CP-even scalars, the MSSM and bilep¬ 
ton sectors are almost decoupled, mixing exclusively due 
to the gauge kinetic mixing. In first approximation, the 
mass matrix is block-diagonal, and has mass eigenstates 
that mimic the MSSM case. In practice, it turns out that 
only two Higgs bosons are light (hereafter called Hi and 
H 2 , one per sector), while the other two are very heavy 
(above the TeV scale). The lightest scalars are well de¬ 
fined states, being either almost exclusively doublet-like 


or bilepton-like. It is worth stressing that their mixing 
is small (see Fig. |4|) and solely due to the gauge kinetic 
mixing (see also Ref. [^b 

Concerning the physical pseudoscalars AP and A®, 
their masses are given by 


m^o — 


2B^ 
sin 2/3 


m^o = 


2B^, 

sin 2/3' 


(17) 


For completeness we note that the mass of charged Higgs 
boson reads as in the MSSM as 


m^+ = Bf^ (tan /3 -|- cot /3) -I-: 


‘-w 


(18) 


In this model, the CP-odd and charged Higgses are 
typically very heavy. In eq. m we see that compared 
to the MSSM, there is a non-negligible contribution from 
the gauge kinetic mixing. LHC searches limit tan/3' < 
1.5 and v-r. > 7 TeV, since 


(19) 


Mz' > 3.5 TeV 


at 95% C.L.. Notice that recent reanalysis of LEP pre¬ 


cision data also constrain fa, > 7 TeV at 99% C.L. [37 1. 


A consequence of this strong constraint in the BLSSM is 
that the first terms in eqs. (USD-iiini) can be large, pushing 
for CP-odd and charged Higgs masses in the TeV range. 

The very large bound on the Z' mass is in contrast 
with the non-SUSY version of the model, where the gauge 
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couplings are free parameters and can be much smaller, 
hence yielding lower mass bounds. The latter need to be 
evaluated as a function of both gauge couplings . 

Next, we describe the sneutrino sector, that shows two 
distinct features compared to the MSSM. Firstly, it gets 
enlarged by the superpartners of the right-handed neutri¬ 
nos. Secondly, even more drastically, a splitting between 
the real and imaginary parts of each sneutrino occurs 
resulting in twelve states: six scalar sneutrinos and six 
pseudoscalar ones 3, 41|. The origin of this splitting is 
the {yf fj term in the superpotential, eq. (ED, which 
is a AL = 2 operator after the breaking of U{1)b-l- In 
the case of complex trilinear couplings or ^-terms, a mix¬ 
ing between the scalar and pseudoscalar particles occurs, 
resulting in 12 mixed states and consequently in a 12 x 12 
mass matrix. 

To gain some feeling for the behaviour of the sneutrino 
masses we can consider a simplified setup: neglecting 
kinetic mixing as well as left-right mixing, the masses of 
the R-sneutrinos at the SUSY scale can be expressed as 

m?s ^ mlc + M|, Q cos( 2/3') -|- sin/3'^^ 

V2Y, 


Mz' 


9bl 


{Ax sin /?' — /r' cos /3') , (20) 


/1 2Y^ 

m~p ~ -I- M|, f - cos(2/3') -|- sin/3'^ 

V 9 BL 

\f^Y 

— Mz' - - {Ax sin P' — fi' cos /3') . (21) 

9bl 


In addition, we treat the parameters Ax, mlc, Mz', 1 *', 
Yx and tan/3' as independent. The different effects on the 
sneutrino masses can easily be understood by inspecting 
eqs. EOD and EID- The first two terms give always a 
positive contribution whereas the third one gives a con¬ 
tribution that can be potentially large which differs in 
sign between the scalar and pseudoscalar states, there¬ 
fore inducing a large mass splitting between the states. 
Further, this contribution can either be positive or nega¬ 
tive depending on the sign of Ax sin /3' — /r' cos /3'. For ex¬ 
ample choosing Yx and fY positive, one finds that the CP- 
even (UP-odd) sneutrino is the lightest one for Ax < Q 
{Ax > 0). This is pictorially shown in Fig. [T1 as a func¬ 
tion of the GUT-scale input parameter Aq, for a choice 
of the other parameters. One notices that the CP-even 
(ClP-odd) sneutrino is the lightest one when the 125 GeV 
Higgs boson is predominantly Hi (i? 2 )- It is worth point¬ 
ing out here that, as will be described in the following sec¬ 
tion, when = 125 GeV, the next-to-lightest Higgs 
boson can decay into pairs of CP-even sneutrinos, but 
not into the similar channel with CP-odd sneutrinos. Be¬ 
ing H 2 predominantly a bilepton field, when this decay 
is open it saturates its BRs, see Fig. [3] Regarding the 
decay into CP-odd sneutrinos, this channel is accessible 
{i.e. is light enough) only in the region where H 2 


is the SM-like Higgs boson, i.e. mainly coming from the 
doublets. In this case however, this decay channel is miti¬ 
gated by the small scalar mixing and is not overwhelming 
(unlike for Hi, now mainly from the bileptons). 



-2000 -1000 0 1000 2000 3000 4000 


A„ [GeV] 

FIG. 1: Masses of CP-even {u^, cyan) and CP-odd {u^, red) 
R-sneutrinos as a function of Aq. For comparison, also the 
masses of the lightest (Pi, black) and next-to-lightest (P 2 , 
blue) Higgs bosons are shown. In green, it is shown configu¬ 
rations when Mhi = 125 GeV. 


Depending on the parameters, either type of sneutri¬ 
nos can get very light. If the LSP, it can be a suitable 
dark matter candidate 16| and yield extra fully invisible 
decay channels to the Higgs bosons, thereby increasing 
their invisible widths. In the case of the decay into the 
CP-odd sneutrino, since this can happen mainly for the 
SM-like Higgs boson, one should account for the con¬ 
straints on the former Eventually, the R-sneutrinos 
could also get tachyonic or develop dangerous P-parity- 
violating vevs. While the first possibility is taken into 
account in our numerical evaluation by SPheno, and such 
points are excluded from our scans, the second case will 
be reviewed in the following subsection. 

The last important sector for considerations that will 
follow is the one of the charged sleptons. See Ref. 
for further details. New SUSY breaking D-term contri¬ 
butions to the masses appear, that can be parametrised 
as a function of the Z' mass and of tan /3' as 


Mz'{tan^ P' - 1) 
2 1 -j- tan^ P' 


( 22 ) 


Their impact is larger for the sleptons than for the 
squarks by a factor of 3 due to the different B — L 
charges {Q^~^). It is possible to vary the stau mass 
by ±0(100) GeV with respect to the MSSM case while 
keeping the impact on the squarks under control. Having 
different sfermion masses in the BLSSM as compared to 
the MSSM has a net impact onto the Higgs phenomenol¬ 
ogy, in particular in enhancing the /177 coupling while 
keeping unaltered the SM-like Higgs coupling to gluons. 
As described at the end of this review, the new D-terms 
coming from the B — L sector can further reduce the 
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stau mass entering in the /177 effective interaction (while 
ensuring a pole mass of ~ 250 GeV, compatible with ex¬ 
clusions) ^ leading this mechanism to work also in the 
constrained version of the model. This mechanism has 
been recently reanalysed also in Ref. 4^ in the very same 
model. 


The issue of R-Parity conservation 


We have encountered so far several neutral scalar fields 
with could develop a vev, beside the Higgs bosons. If 
vevs of fields charged under QCD and electromagnetism 
are forbidden because the latter are good symmetries, 
R-sneutrino vevs, which are not by themselves problem¬ 
atic, would unavoidably break R-Parity. The issue of 
conserving R-Parity is of fundamental importance, since 
this is a built-in symmetry in our model where B — L is 
gauged. We will therefore restrain ourselves to parame¬ 
ter configurations where the global minimum is R-Parity 
conserving. 

When all neutral scalar fields are allowed to get a vev, 
it is not trivial even at the tree level to find which is the 
deeper global minimum, and whether it is of a “good” 
type, here defined as having the correct broken symme¬ 
tries and being R-Parity conserving. One possible way to 
study this issue is to start from a simplified set of input 
parameters yielding a correct tree level global minimum 
when only the Higgs fields get a vev, and then look for 
the true global minimum when all other neutral fields 
(mainly R-sneutrinos) acquire a vev, both at the tree 
level and at loop level. See Ref. ^ for further details. 

At the tree level there seems to exist regions where the 
BLSSM has a stable, R-Parity-conserving global mini¬ 
mum with the correct broken and unbroken gauge groups. 
For this to happen one needs the R-sneutrino Yukawa 
coupling Yx to be not so large, and the trilinear param¬ 
eter Aq to be not large compared to the soft scalar mass 
mo, as, intuitively, large and Aq can lead to large 
negative contributions to the potential energy for large 
values of Vx , as well as reducing the effective R-sneutrino 
masses, as described above and clear from Fig. [TJ 

It turns out that when loop corrections are taken into 
account, few points all over such regions of parame¬ 
ters exist where R-Parity is not preserved anymore, or 
where SU{2)l or U{1)b-l are unbroken. This is appar¬ 
ently due to a very finely-tuned breaking of SU{2)l and 
U{1)b-l which often does not survive loop corrections. 
The reason for this is that besides the known large con- 


^ With pole mass we denote the one-loop corrected mass at Q = 
^SUSY = \/tit 2 , while in the loop, leading to the effective h-yy 
coupling, the running DR tree-level mass at Q = enters, 
being h the SM-like Higgs boson, i.e. mj, = 125 GeV. 


tributions of third generation (s)fermions, the additional 
new particles of the B — L sector also play an impor¬ 
tant role. As previously for the charged sleptons sector, 
new SUSY breaking D-term contributions to the masses 
appear, see eq. (1^^ . Since, as shown in eq. (USD, the ex¬ 
perimental bounds require Mz' to be in the multi-TeV 
range, these contributions can be much larger than in 
the MSSM sector, resulting in the observed importance 
of the corresponding loop contributions. Furthermore, 
these contributions are also responsible for the restora¬ 
tion of U{1)b-l at the one-loop level. 

Ultimately, overall safe regions of parameters cannot 
be found where the correct vacuum structure can be en¬ 
sured. At the same time, if naive trends can be spotted 
for bad points to appear, these have nonetheless to be 
checked case-by-case due to the highly non-trivial scalar 
potential, and it might be possible that neighbour config¬ 
urations still hold a valid global minimum. We will not 
check the validity of our scans from the vacuum point of 
view in the following, being confident that if any point 
is ruled out, a neighbour one yielding a very similar phe¬ 
nomenology can be found, which is allowed. 


A QUICK LOOK TO FLAVOUR OBSERVABLES 


Before moving to the Higgs phenomenology, we briefly 
show the impact on the BLSSM model when considering 
the constraints arising from low energy observables. For 
a review of the observables as well as for the impact onto 
general SUSY models encompassing a seesaw mechanism, 
see Refs. 45, 46|. 


We consider here only the two most constraining ones, 
BR(/r —>• 67 ) and BR(/r —>• 3e). The present exclusions 
are BR(^ —f 07 ) < 5.7 • 10“^^ [i^ and BR(/r —f 3e) < 
1 • 10“^^ [ 4 ^. In Fig. [2] we plot these branching ratios 
as a function of the mass of the lightest (in black) and 
next-to-lightest (in red) SM-like neutrino, which display 
some pattern for evading the bounds. In particular, they 
are required to be rather light, below 0.5 eV, while the 
model, ought to the scans here performed, seems to prefer 
configurations with neutrinos heavier than O.OI eV, hence 
the preferred region in between. Lighter mass values are 
nonetheless also allowed. 


For convenience, the impact of satisfying the earlier 
bounds will be shown only in the inverted hierarchy case, 
due to the smaller density of configurations therein. In¬ 
stead, points not allowed in the normal hierarchy case 
are automatically dropped. 

Regarding the long-lasting {g — 2)^ discrepancy, in the 
setup investigated here charginos and charged Higgses are 
too heavy, same for the Z' boson, while the neutralino 
and sneutrino are too weakly coupled, to give a significant 
enhancement over the SM prediction. 
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due to the kinetic mixing. 



FIG. 3: Branching ratios for H 2 with Mh 2 > Mhi = 125 
GeV. The CP-even sneutrino channel (brown) is superim¬ 
posed. 


FIG. 2: (Upper plot) BR(/r —>■ 67 ) and (lower plot) BR(/r —>■ 
3e) as a function of the light neutrino masses in GeV (black: 
vi, red: 122 )■ The blue horizontal lines represent the ac¬ 
tual experimental limits, from Refs, and [i^ . respec¬ 
tively. The parameters have been chosen as mo € [0.4, 2] TeV, 
Mi /2 e [1.0,2.0] TeV, tan/3 G [5,40], £ [-4.0,4.0] TeV, 

tan/3' G [1.05,1.15], M^/ G [2.5, 3.5] TeV, ^ G 1 • [0.002, 0.4], 
n G 1 ■ [0.05,5] X 10"®. 


HIGGS PHENOMENOLOGY 


We review here the phenomenology of the Higgs sec¬ 
tor, showing a first survey of its phenomenological fea¬ 
tures. First, results when normal hierarchy is imposed 
are presented. Then, we will show that the inverted hier¬ 
archy is also possible on a large portion of the parameter 
space. Without aim for completeness, the results are here 
presented as the starting point for a more thorough in¬ 
vestigation. Finally, it is described how model features 
pertaining to the extended gauge sector impinge onto the 
Higgs phenomenology, and in particular how the Higgs- 
to-diphoton branching ratio can be easily enhanced in 
this model, despite the experimental data now converg¬ 
ing to a more SM-like behaviour than in the recent past. 


Normal hierarchy 

In this subsection we discuss the normal hierarchy case, 
with the lightest Higgs boson being the SM-like one {i.e., 
predominantly from the doublets), and a heavier Higgs 
boson predominantly from the bilepton fields (those car¬ 
rying B — L number and responsible for its spontaneous 
breaking). Their mixing is going to be small and solely 


In Fig. Owe first inspect the heavy Higgs boson branch¬ 
ing ratios. Besides the standard decay modes, the decay 
into a pair of SM Higgs bosons exist, as well as two new 
characteristic channels of this model, comprising right- 
handed (s)neutrinos. 

1. H 2 —>■ HiHi. Its HR can be up to 40% before the 
top quark threshold, and around 30% afterwards; 

2 . H 2 VhVh- A similar decay channel exists for 
the Z' boson. The HR are 0(10)%, up to 20% 
depending on the heavy Higgs and neutrino masses; 

3. H 2 —>■ where, is the CP-even sneutrino 

and the LSP, hence providing fully invisible decays 
of the heavy Higgs. If kinematically open, it satu¬ 
rates the Higgs BRs. Notice that only points with 
very light CP-even sneutrinos are shown, possible 
only for very large and negative Aq (see Fig. [T|). 


While the first two channels exist also in the non-SUSY 
version of the model ^ (see, e.g., 0), the last one, in¬ 
volving the CP-even sneutrino, is truly new and rather 
intriguing. This is because the sneutrino is light and it 
can be a viable LSP candidate if with mass lower than 
H 2 , as in this case . It however implies that the heavy 
Higgs is predominantly bilepton-like, with a light Higgs 
very much SM-like. This can be seen in Fig.|4l where the 
points with large BR(iL 2 —>■ (in red) have the low¬ 

est mixing between H 2 and the SM scalar doublet fields, 
of the order of 0.1%. It immediately follows that this 


® However, in the non-SUSY B — L model the Higgs mixing angle 
is a free parameter, directly impacting on these branching ratios. 
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150 200 250 300 350 400 450 

Ma, [GeV] 


FIG. 4: Mixing between Higgs boson mass eigenstates (blue- 
orange: Hi, cyan-red: H 2 ) and scalar doublet fields, as a 
function of Mh 2 - ZH[i,j] is the scalar mixing matrix. Or¬ 
ange/red points are the subset corresponding to BR(772 —>■ 
> 90% . 


channel will have very small cross section at the LHC, 
when considering SM-like Higgs production mechanisms. 
This is true for all heavy Higgs masses Mh 2 > 140 GeV. 
The 125 GeV Higgs is well SM-like, with tiny reduction 
of its couplings to the SM particle content. On the other 
side, the heavy Higgs is feebly mixed with the doublets, 
suppressing its interactions with the SM particles, and 
hence its production cross section. This can be seen in 
Fig. [ 5 ] (top frame). Considering only the gluon fusion 
production mechanism, and multiplying it by the rele¬ 
vant HR, we get the cross sections for the choice of chan¬ 
nels displayed therein. The most constraining channels, 
H —>• WW —f (I'jj and H —>• WW —f 2£2i', are also com¬ 
pared to the exclusions at the LHC for ^ = 8 TeV from 
Refs. [53 and M|, respectively. The H —>• ZZ chan¬ 
nels are well below current exclusions, that are hence not 
shown. 


We see that all the displayed configurations are al¬ 
lowed by the current searches (the exclusions shown by 
solid curves of same color as the depicted channel). This 
is because of the suppression of the heavy Higgs boson 
cross sections due to the small scalar mixing. 

In the lower plot are displayed the cross sections for 
the new channels. Those pertaining to model configu¬ 
rations for which the heavy Higgs boson decays to the 
CF-even sneutrino (LSP), yielding a fully invisible decay 
mode, are displayed in red. Contrary to the all other 
cases, the production of the heavy Higgs for this channel 
is via vector boson fusion as searched for at the LHC . 
Typical cross sections range between 0.1 fb and 1 fb. The 


Starting from > 130 GeV. 



FIG. 5: Gross sections at ^/s = 8 TeV for (upper plot) the SM- 
like channels (lower plot) the new channels, as a function of 
the heavy Higgs mass. The solid lines above are the exclusion 
curves from [Hfil . [Ml- 


H 2 —f HiHi channel is shown in blue and it can yield 
cross sections of 1 4- 10 fb for 250 < Mh 2 < 400 GeV. 
Last is the H 2 —f Vh^h channel. It can be sizable only for 
very light H 2 masses: ~ lO-FlOO fb for 140 < Mh 2 < 160 
GeV, although the further decay chain of the heavy neu¬ 
trinos have to be accounted for. The latter can give spec¬ 
tacular multi-leptonic final states of the heavy Higgs bo¬ 
son {A£2v and M2jv) or high jet multiplicity ones (2I'4j), 
via Vh —f and Uh —f vZ in a 2 : 1 ratio (mod¬ 

ulo threshold effects). Further, these decays are typically 
seesaw-suppressed and can therefore give rise to displaced 
vertices [^. 


Inverted hierarchy 

In this subsection we discuss the inverted hierarchy 
case, where H 2 is the SM-like boson and a lighter Higgs 
boson exists. 

We start once again by presenting the HRs for the 
next-to-lightest Higgs boson in Fig. [51 This time how¬ 
ever this is the SM-like boson, hence predominantly from 
















in this case displaced vertices may appear. 


0.1 .TT 


• bb ; 

■iv±vr=F- 





125.0 125.5 

Mh, [GeV] 


FIG. 6: Branching ratios for the 125 GeV Higgs boson ( 1 ^ 2 ). 
The decay into heavy neutrinos is displayed with diamonds. 
All others with circles. Gray points are excluded by the low 
energy observables and by HiggsBounds. The decay into CP- 
odd sneutrinos is not shown. 


the doublets. It has the same new channels as the heavy 
Higgs in the normal hierarchy, the only difference be¬ 
ing the CF-odd R-sneutrino instead of the CF-even one. 
This is simply because the inverted hierarchy can hap¬ 
pen only for large positive Aq values, where only the 
CF-odd R-sneutrino can be light, see Fig. [T] The con¬ 
figurations not allowed by the low energy observables or 
by HiggsBounds are displayed as gray points. We see 
that H 2 may have sizable decays into pairs of the lighter 
Higgs bosons, yielding 4Fjets final states. This decay is 
still allowed with rates up to few percent. Further, rare 
decays into pairs of heavy neutrinos are also present, with 
BRs below the permil level. This channel can give rise 
to rare multi-lepton/jets decays for the SM-like Higgs 
boson, that are searched for at the LHC, even in com¬ 
bination with searches for displaced vertices 5^. The 
last available channel is the decay into pairs of CF-odd 
R-sneutrinos. Being the LSP, it will increase the invisible 
decay width and hence give larger-than-expected widths 
for the SM-like boson. Its rate is obviously constrained, 
and a precise evaluation of the allowed range is needed. 
It however goes beyond the scope of the present review 
and we postpone it to a future publication. 


Regarding the lightest Higgs boson (Hi), this will ob¬ 
viously decay predominantly into pairs of 6-jets. Notice 
that due to its large bilepton fraction it can also decay 
into pairs of very light RH neutrinos, at sizable rates 
depending on the neutrino masses. As in the in previ¬ 
ous figure, the non-allowed configurations are displayed 
as gray points. We see that the pattern of decays is not 
affected by the inclusion of the constraints, in the sense 
that this channel stays viable. Once again, the latter will 
yield multi-lepton/jet final state, which will be very soft, 
and hence very challenging for the LHC. However, also 



FIG. 7: Same as in Fig. [6]for the lightest Higgs boson (Hi). 


As in the previous section, we show in Fig. [5] the mix¬ 
ing between the Higgs mass eigenstates and the doublet 
fields as a function of the light Higgs mass, to show that 
H 2 is here rather SM-like. Once more, the gray points 
displayed here are excluded by the low energy observables 
and by HiggsBounds. 



20 40 60 80 100 120 

[GeV] 


FIG. 8: Mixing between scalar mass eigenstates and Higgs 
doublets, (black: H\, red: H 2 ) and scalar doublet fields, as a 
function of Mhi- ZH[i,j] is the scalar mixing matrix. Gray 
points are excluded by the low energy observables and by 
HiggsBounds. 


Finally, the production cross sections for the lightest 
Higgs boson can be evaluated. In Fig. Owe compare the 
direct production (for the main SM production mecha¬ 
nisms, gluon fusion and vector boson fusion) with the 
pair production via H 2 decays only via gluon fusion, 
gg H 2 —?► Hi Hi. When the latter channel is kine¬ 
matically open, i.e. 2 Mh^ < 125 GeV, the lightest Higgs 
boson pair production has cross sections up to 1 pb at 











9 


the LHC at y/s = 8 TeV, and it can give rare 46, 2b2V 
or 4V (V = W, Z) decays of the SM-like Higgs boson. 
A thorough analysis of the phenomenology of the Higgs 
sector in the BLSSM for the upcoming LHC run 2, based 
on the first investigations shown here, will be performed 
soon. 




b 
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0.001 


20 40 60 80 100 120 


Mh, [GeV] 


corresponding to the contributions from charged SM 
fermions, W bosons, charged Higgs, charginos, charged 
sleptons and squarks, respectively. The amplitudes Ai 
at lowest order for the spin-1, spin-i and spin-0 parti¬ 
cle contributions, can be found for instance in Ref. [55| . 
ghxx denotes the coupling between the Higgs boson and 
the particle in the loop and Qx is its electric charge. In 
the SM, the largest contribution is given by the IH-loop, 
while the top-loop leads to a small reduction of the decay 
rate. In the MSSM, it is possible to get large contribu¬ 
tions due to sleptons and squarks, although it is difficult 
to realise such a scenario in a constrained model with 
universal sfermion masses In singlet or triplet 

extension of the MSSM also the chargino and charged 
Higgs can enhance the loop significantly [^, . How¬ 

ever, this is only possible for large singlet couplings which 
lead to a cut-off well below the GUT scale. In contrast, it 
is possible to enhance the diphoton ratio in the BLSSM 
due to light staus even in the case of universal boundary 
conditions at the GUT scale. We show this by calculating 
explicitly the contributions of the stau: 


FIG. 9: Cross sections at ^/s = 8 TeV for different production 
mechanisms. Gluon-fusion (in red) and vector-boson-fusion 
(in green) mechanisms are displayed only for Mhi > 50 GeV 
for simplicity. Gray points are excluded by the low energy 
observables and by HiggsBounds. 


1 9detm| 

3 cHogn 

2 2m^(AT-—/itan/3)^ 

3 (m|; -I- D}i){m\ -b Dl) + m’^gtanP{2Ar 


(24) 


p,tan P) 
(25) 


(see, e.g., [^) 


- h—¥'yy — 




128^271 


Here, and Du represent the D-term contributions of 
the left- and right-handed stau and we have neglected 
sub-leading contributions. Given that 2Ar < /itan,5, for 
fixed values of the other parameters, Du and D^ can be 
used to enhance the 77 rate by suppressing the denomi¬ 
nator. 

We turn now to a fully numerical analysis to demon¬ 
strate the mechanism to enhance the Higgs-to-diphoton 
rate as a feature of the model with an extended gauge 
sector. This is a result of reducing the stau mass at 
the Higgs mass scale via extra D-terms as shown dis¬ 
cussing eq. (I22I1 . We remind here that this mechanism 
leaves the stop mass and hence, as we will show, the 
Higgs-to-gluons effective coupling nearly unchanged. In 
Table m we have collected two possible scenarios that 
provide a SM-like Higgs particle in the mass range pre¬ 
ferred by LHG results displaying an enhanced diphoton 
rate. In the first point, the lightest CP-even scalar eigen¬ 
state is the SM-like Higgs boson while the light bilepton 
is roughly twice as heavy. In Fig. [in] we show that all 

E 2 Ah ( \ Ah( features arise from the extended gauge sector: it is 

cQfQhff i/ 2 {'''f) + 9hww 1 G^lgufficient to change only tan/3' to obtain an enhanced 

77)]b-i 


Enhancement of the diphoton rate 

A feature of gauge-extended models is that new SUSY- 
breaking D-terms arise, that give further contributions 
to the sparticle masses. In the case of the model under 
consideration, we showed discussing eq. (1221) that these 
terms can be large, and that they bring larger corrections 
to sleptons than to squarks. We already discussed how 
the vacuum structure of the BLSSM is affected by this. 
Here, we discuss their impact on the Higgs phenomenol¬ 
ogy, focusing on the Higgs-to-diphoton decay, despite dis¬ 
favoured W most recent data [17|, as an illustrative case. 
See Ref. for further details. 

To start our discussion let us briefly review the partial 
decay width of the Higgs boson h into two photons within 
the MSSM and its singlet extensions. This can be written 
as 


/ 


-b 


hny^ghH+H- 








E 


2ray/ 


diphoton signal 


^xt 


Ah 

f^X+X” 1/2 


E QhejCj 

?T7? 




E 


ghqiq 


■iQlKirq.) 




(23) 


and the 


SM 

correct dark matter relic density while keeping the mass 
of the SM-like Higgs nearly unchanged. The dark matter 
candidate in this scenario is the lightest neutralino, that 
is mostly a bileptino (the superpartner of the bileptons). 
The correct abundance for tan /3' ~ 1.156 is obtained due 
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1.140 1.145 1.150 1.155 


tan /?' 

FIG. 10; (Top plot) The mass of the SM-like Higgs [bot- 
tom(blue line)], of the stau [middle(black) line, where the 
dashed line represents a reference unchanged value] and of the 
lightest neutralino [top(red) line]; (middle plot) the diphoton 
branching ratio; (bottom plot) the neutralino relic density as 
a function of tan/3'. The other parameters have been cho¬ 
sen as mo = 673 GeV, M- 1/2 = 2220 GeV, tan/3 = 42.2, 
^0 = -1842.6, Mz' = 2550 GeV, K, = 1 • 0.42 


to a co-annihilation with the light stau. In the second 
point, the SM-like Higgs is accompanied by a light scalar 
around 98 GeV which couples weakly to the SM gauge 
bosons, compatibly with the LEP excess (bll - l^ . In this 
case, the LSP is a CP-odd sneutrino which annihilates 
very efficiently due to the large Y^- This usually results 
in a small relic density. To get an abundance which is 
large enough to explain the dark matter relic, the mass 
of the sneutrino has to be tuned below mw [l^ . This 
can be achieved by slightly increasing tan P' and by tun¬ 
ing the Majorana Yukawa couplings that tends to 
increase the SM-like Higgs mass for the given point. It is 
worth mentioning that a neutralino LSP with the correct 
relic density in the stau co-annihilation region can also be 
found in this scenario. Notice that both points yield rates 
consistent with observations in the WW* jZZ* channels 
(measured at the LHC) (being Chzz ^ 1), as well as an 
effective Higgs-to-gluon coupling close to I. 



Point I 

Point 11 

mjtj [GeV] 

125.2 

98.2 

mh 2 [GeV] 

186.9 

123.0 

ruf [GeV] 

267.0 

237.3 

doublet fr. [%] 

99.5 

8.7 

bilepton fr. [%] 

0.5 

91.3 


0.992 

0.087 

ChiZZ 

1.001 

0.085 

^h299 

0.005 

0.911 

Ch2ZZ 

0.005 

0.921 

r(/ii) ]MeV] 

4.13 

0.22 


1.57 

0.085 

Rh 

1.03 

0.089 

Rww* 

0.98 

0.05 

r(/i2) [MeV] 

4.8 

3.58 


0.005 

1.79 

Ki 

0.006 

0.95 

Rww” 

0.01 

0.88 

LSP mass ]GeV] 

253.9 

82.9 

Q.h^ 

0.10 

10"^ 


TABLE II: The input parameter used: Point I: mo = 
673 GeV , M 1/2 = 2220 GeV, Ao = -1842 GeV, tan ,3 = 42.2, 
tan/3' = 1.1556, = 2550 GeV, Y, = 1 ■ 0.42 (neutralino 

LSP). Point II: mo = 742 GeV , M 1/2 = 1572 GeV, Ao = 
3277 GeV, tan/3 = 37.8, tan ,3' = 1.140, M^/ = 2365 GeV, 
Yx = diag(0.40, 0.40, 0.13) (CP-odd sneutrino LSP). csvv de¬ 
notes the coupling squared of the Higgs fields to vector bosons 
normalised to the SM values. 

CONCLUSIONS 

In this review I described the U{\)b-l extension of 
the MSSM, focusing in particular on the scalar sector, 
described in details. The fundamental role that the gauge 
kinetic mixing plays in this sector has been underlined. 

The comparison to the most constraining low energy 
observables showed that a preferred region for the light 
neutrino masses exists to evade these bounds. Then, 

1 presented a first systematic investigation of the phe¬ 
nomenology of the Higgs sector of this model, showing 
that both the normal hierarchy and the inverted hierar¬ 
chy of the two lightest Higgs bosons are naturally pos¬ 
sible in a large portion of the parameter space. Par¬ 
ticular attention has been devoted to analyse the new 
decay channels comprising both the CP-even and CP- 
odd R-sneutrinos, which are a peculiarity of the BLSSM. 
Based on these first findings, a thorough analysis of the 
Higgs sector in the BLSSM at the upcoming LHC run 

2 will be soon prepared. The ht of the SM-like Higgs 
boson to the LHC data will also be performed with 
HiggsSignals [b^. 

Finally, I described how in the BLSSM model (and 
in general in gauge-extended MSSM models) the Higgs- 
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to-diphoton decay can be easily enhanced. Despite dis¬ 
favoured by most recent data, this feature is a conse¬ 
quence of the potentially large new SUSY-breaking D- 
terms arising from the B — L sector. At the same time 
these terms affect also the vacuum structure of the model, 
where naive R-Parity conserving configurations at the 
tree level, could develop deeper R-Parity violating global 
minima, or partially restore the SU{2)l x U{1)b-l sym¬ 
metry at one loop. It is however possible to still find R- 
Parity conserving global minima on the whole parameter 
space, which can either accommodate an enhancement of 
the Higgs-to-diphoton decay or fit the most recent Higgs 
data. 
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